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We study the scaling properties in deep inelastic scattering using the most recent combined struc- 
ture function data F2 from the HI and ZEUS collaborations. We also perform a direct fit to the 
F2 data inspired by the scaling properties. Our analysis favours the QCD saturation mechanism 
from the Balitski Kovchegiv equation wuth running coupling. 
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1. Scalings in deep inelastic scattering 

Geometric scaling [p], §1 is a remarkable empirical property found originally using the data 
on high energy deep inelastic scattering (DIS) i.e. virtual photon-proton cross-sections. One can 
represent with reasonable accuracy the cross section a r * p by the formula a' y * p (Y, Q) = (t) , 
where Q is the virtuality of the photon, Y the total rapidity in the /"-proton system and 

T = io g <2 2 -iog<2 J (F)=iog<2 2 -Ar , (i.i) 

is the scaling variable. A fit to the DIS data measured leads to a value of A ~ .3, which confirms 
the value found within the Golec-Biernat and Wusfhoff model [||] where geometric scaling was 
explicitely used for the parametrization. 

The scaling using the variable x defined in Formula [O] is directly related to the concept of 
saturation, the behaviour of perturbative QCD amplitudes when the density of partons becomes 
high enough. There were many theoretical arguments to infer that in a domain in Y and Q 2 where 
saturation effects set in, geometric scaling is expected to occur. Within this framework, the function 
Q S {Y) can be called the saturation scale, since it determines the approximate upper bound of the 
saturation domain. 

This type of geometric scaling is motivated by asymptotic properties of QCD evolution equa- 
tions with rapidity. Using the nonlinear Balitsky-Kovchegov (BK) equation [Q] which represents 
the "mean-field" approximation of high energy (or high density) QCD, geometric scaling could 
be derived from its asymptotic solutions j^j]. This equation is supposed to capture some essential 
features of saturation effects. Considering the BK equation with fixed coupling constant leads 



asymptotically to the original geometric scaling of Formula 1.1. Considering a running coupling 
leads to the following scaling 

T = log Q 2 - log Q s (Y ) = log Q 2 - X VY , (1.2) 



Recently flBh, it was noticed that the scaling solution [L2| of the BK equation with running 
coupling is only approximate and not unique. Another equivalent approximation leads to a different 
scaling variable, namely 

T = logG 2 -A — . (1.3) 

The effect of QCD fluctuations was examined in Ref. [j7|] in the fixed coupling scheme and 
gives rise to a new "diffusive scaling", the scaling variable being 

log2 2 -AF 

T= SV V • (1.4) 

The aim of this study from a theoretical point-of-view is to test and compare the different 
scaling behaviors, arising from different versions of QCD evolution, using the most recent precise 
data available from HERA resulting from a combination of the HI and ZEUS F2 measurements [j|]. 
We study the quality of the description of these combined data set using the four kinds of scaling 



and refer to them in the following as "Fixed Coupling" for the variable 1.1, "Running Coupling I" 
and "Running Coupling II" for the variables [L2| and [O] respectively, and to "Diffusive Scaling" 
forO. 
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Running Coupling I 
O HERA Data 



Figure 1: DIS cross section data as a function of the scaling variable for Running Coupling I. 



2. The Quality Factor 

In order to compare the quality of the different scalings and to check if the DIS cross sections 
~ F2/Q 2 depend mainly on the T variable or not, it is useful to introduce the concept of Quality 
Factor [^] (QF) while the explicit form of the T dependence is not known. After normalising the 
data sets v,- = log (a, ) and scalings Uj = T, (A) between and 1, and ordering the scalings in it,-, we 
define QF 



QF{X) 



(Vi~Vi-l) 2 
: ( --W;_i) 2 + £ 2 



e is needed in the case that two data points have the same scaling, namely when they have the same 
x and Q 2 , and we take £=0.01. The method is to fit the value of A to maximize QF. 



3. Scaling tests in DIS 

As we mentioned, we use the very precise data sets combining the HI and ZEUS measure- 
ments of the proton structure function Fi [Js[] - To remain in the region where perturbative QCD is 
applicable and to avoid the region where valence quarks dominate, we choose to restrict ourselves 
to data points with 4 < Q 2 < 150 GeV 2 and x < 10~ 2 . In addition, in order to avoid the high y 
region where Fl is large, we add an additional cut on data on y < 0.6. After all cuts, we are left 
with 117 data points. 

The values of the X parameters and the QF are given in Table I for the different scaling consid- 
ered in this analysis. While Fixed Coupling, Running Coupling I and II lead to approximately the 
sames value of QF, Diffusive Scaling is clearly disfavoured. As an example, the scaling plot show- 
ing all combined DIS cross section data as a function of the X variable for Running Coupling I is 
given in Fig. I to show the quality of scaling. In addition, it is worth noticing that adding additional 
variables such as Qq or a shift in rapidity Yq does not improve the scaling quality. 
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scaling 


parameter 


1/QF 


Fixed Coupling 


A = 0.31 


150.2 


Running Coupling I 


A = 1.61 


137.9 


Running Coupling II 


A =2.76 


124.3 


Diffusive Scaling 


A =0.31 


210.7 



Table 1: Values of 1 / QF and of the A parameter for the four different kinds of scaling considered 




Figure 2: Result of the Running Coupling I fit 
to the combined Fo data set from HI and ZEUS. 
We note a fair description of data for x < 10~ 2 
and y < 0.6. 



Figure 3: Extrapolation of the Running Cou- 
pling I fit to the combined Fi data set from HI 
and ZEUS at low Q 2 values. 



4. Fits to HERA data 



In this section, we describe a fit to the combined HERA data motivated by the success of the 
data description using the Running Coupling I scaling variable. In the fit, we will use all data 
above Q 2 = AGeV 2 since the fitting formula that we develop is valid only in the dilute regime, and 
saturation is supposed to occur at very low Q 2 at HERA. The following formula, deduced from the 
dipole amplitude with saturation including the asymptotic expression of the Airy function which is 
the solution of the Balistsky-Kovchegov equation, is used to fit the data 

T = log&)-A A /log(i)-y 

( -/3t 3 /2 \ 
a = iVexp(— ar)exp — - 

yy ' y \(iogi/x-Y y/*) 

where the different parameters used in the fit are A, a, j8, Qo, Yq and N. We notice that this formula 
shows only a moderate scaling violation introduced by the (log l/x — Yq) 1 / 4 term predicted by the 
dipole model and we perform the fits with and without this term. 

The fit results and the parameter values are given in Table II and Fig. II. We note that the fit 
X 2 is close to 1.2 per dof and is similar with or without the scaling violation term. The fit cannot 
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Parameter 


Fit I 


Fit II 


X 


1.54 ±0.02 


1.54 ±0.02 


a 


0.34 ± 0.01 


0.18 ±0.01 


P 


0.24 ± 0.01 


0.18 ±0.01 


Go 


0.079 ± 0.01 


0.064 ± 0.01 


Y 


-1.46 ±0.02 


0.50 ± 0.02 


N 


0.51 ±0,01 


0.72 ± 0,01 


X 2 


130.1 


119.0 



Table 2: Value of the parameters of the Running Coupling I inspired fit to the combined DIS cross section 
1 17 data points from HERA. Fit I (resp. Fit II) is performed with (resp. without) the scaling violation term. 

describe the reduction of the reduced cross section at high y due to the large values of Fl. In Fig. Ill, 
we also show the fit extrapolation at lower Q 2 which leads to a fair description of data. Going to 
lower values of Q 2 will require a parametrisation valid in the saturated region whereas our formula 
is only valid in the dilute regime. 

In addition, we also attempted to perform a similar fit inspired by Fixed Coupling or Running 
Coupling II, but they lead to a worse description of data (j 2 =156.4 and 190.4 respectively). 

5. Conclusion 

In this paper, we presented a new study of scaling properties in DIS using the most recent 
combined F2 data from HI and ZEUS. The new precise data set shows precise scaling using either 
Fixed Coupling, Running Coupling I or II, while Diffusive Scaling is disfavoured. A direct fit 
inspired by Running Coupling I leads to a good description of data in the dilute regime. 
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